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Subj:  Workshop  on  Acoustic-Propagation  Modeling  by  Non-Ray 
Tracing  Techniques 

1.  The  Acoustic  Environmental  Support  Detachment  (AESD)  was 
established  in  July  1972  as  an  Advanced  Development  activity  to 
synthesize,  evaluate,  document,  and  distribute  a  family  of  Navy 
standard  acoustic  models.  One  of  the  methods  utilized  by  AESD 
to  keep  abreast  of  promising  new  techniques  in  the  Research  and 
Exploratory  Development  communities  is  the  sponsorship  of  workshops 
to  address  particular  aspects  of  the  acoustic  modeling  process. 

2.  This  Technical  Note  provides  a  synopsis  of  such  a  workshop 
which  addressed  the  prediction  of  transmission  loss  by  techniques 
which  did  not  involve  ray-tracing.  It  Is  being  distributed  to 
disseminate  to  a  much  broader  audience  a  summary  of  the  ideas 
presented. 

3.  The  workshop  owes  its  success  to  the  enthusiastic  participation 
of  those  who  attended.  Their  contribution  and  cooperation  are 
gratefully  acknowledged. 


P.  R.  TATRO 
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A  SYNOPSIS  OF  THE  AESD  WORKSHOP  ON 
ACOUSTIC-PROPAGATION  MODELING  BY  NON-RAY-TRACING  TECHNIQUES 
22-25  MAY  1973,  WASHINGTON,  D.  C. 


INTRODUCTION 


From  22  to  25  May  1973  the  Acoustic  Environmental  Support  Detachment 
(AESD)  of  the  Office  of  Naval  Research  sponsored  a  workshop  on  non-ray¬ 
tracing  techniques  used  to  model  acoustic  propagation  in  the  ocean 
environment.  The  basic  objectives  of  the  workshop  were  twofold: 

1.  To  bring  together  those  in  the  R&D  community  working  on  non¬ 
ray  approaches  for  an  informal,  though  critical  review  of  the  current 
state-of-the-art. 

2.  To  indicate  to  the  R&D  community  the  requirements  of  the 
Navy's  Advanced  Development  Program  (as  implemented  by  AESD)  in  this 
field,  thereby  illuminating  new  areas  of  concentration  and  research  to 
meet  future  Navy  needs. 

In  order  to  maintain  a  workshop  atmosphere  and  conform  to  space 
limitations,  attendance  was  limited  to  those  with  working  models  who 
were  willing  to  exercise  them  against  one  or  more  test  cases  provided 
by  AESD.  An  additional  small  group  of  researchers  with  recognized 
credentials  in  the  field  was  invited  to  participate  in  an  advisory 
capacity.  The  workshop  was  divided  into  two  basic  sections:  the  pre¬ 
sentation  by  each  participant  of  his  model  (with  some  discussion  of 
the  predictions  for  each  test  case)  and  workshop  discussions  of  a 
number  of  topics  of  current  interest  relating  to  model  applications. 

It  is  the  purpose  of  this  note  to  report  on  these  proceedings  for 
the  benefit  of  both  the  participants  and  those  who  were  unable  to  attend. 
This  note  consists  of  five  basic  sections: 

1.  A  brief  description  of  the  test  cases. 

2.  A  list  of  attendees  (participants  and  advisors),  the  models 
presented,  and  the  results  for  the  test  cases  against  which  each  mode!- 
was  exercised. 

3.  A  background  discussion  of  solutions  to  the  wave  equation 
to  serve  as  a  basis  for: 

4.  A  brief  description  of  the  approaches  contained  in  each 
model  with  emphasis  on  the  differences  between  models. 

5.  Summaries  of  working  sessions  on: 

a.  Rough-Surface  Scattering 

b.  Volume  Attenuation 

c.  Range  Dependence  of  the  Environment 
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d.  Modeling  the  Ocean  Bottom 

e.  Modeling  Source  Characteristics 

f.  Smoothing  Wave  Results 

g.  Signal-Processing  Applications 


TEST  CASES  AND  PARTICIPANTS 


Three  basic  test  cases  were  provided  to  determine  the  various  regimes 
of  applicability  of  each  model.  In  each  case  the  medium  is  modeled  as 
stratified  (no  range  dependence)  with  a  simple  two-fluid  condition  at 
the  ocean  bottom  characterized  by  a  pc  or  impedance  discontinuity. 

Figure  1  illustrates  Test  Case  1:  characterization  of  a  typical 

North  Pacific  sound  speed  profile  adjusted  to  include  a  very  deep  sur¬ 
face  duct.  Case  1A  is  for  a  source  and  receiver  in  the  duct  at  a  fre¬ 
quency  of  300  Hz  where  two  modes  are  expected  to  be  trapped.  Case  IB 
corresponds  to  a  deep-water  long-range  case.  The  plane-wave  intensity 
reflection  coefficient  for  the  bottom  as  a  function  of  grazing  angle  is 
plotted  in  the  lower  right.  The  numerical  values  for  this  case  are 
shown  in  Figure  2. 

Test  Case  2  is  described  in  Figures  3  and  4  and  is  characteristic 
of  the  double  sound-channel  structure  found  in  the  Northeast  Atlantic. 

Test  Case  3  (Figures  5  and  6)  is  a  shallow  water  case  with  a  relatively 
high-loss  bottom.  At  500  Hz  several  modes  are  trapped,  while  at  50  Hz 
only  1  mode  is  trapped.  The  low-frequency  case  was  chosen  to  test  the 
versatility  of  the  models,  and,  as  the  subsequent  comparisons  will  show, 
the  differences  between  model  predictions  for  this  case  are  the  largest. 

Predictions  generated  by  16  different  models  were  received  for  one 
or  more  of  the  test  cases.  Figure  7  lists  the  models,  the  participants 
or  principal  investigators  responsible  for  the  models,  and  the  cases 
for  which  results  were  received.  The  models  were  divided  into  four 
basic  categories  and  were  presented  in  the  order  shown.  The  parentheses 
around  the  APL  and  CU  cases  indicate  that  normalization  difficulties 
made  comparisons  impossible  and  they  have  not  been  included  in  this 
synopsis.  The  NUSC-3  result  for  test  case  3B  was  off  scale  and  also  has 
been  omitted.  H.  Kutschale  of  Lamont-Dougherty  did  not  receive  an 
invitation  in  time  to  prepare  the  test  cases;  F.  DiNapoli  made  the 
model  presentation  for  W.  Kanabis,  who  was  unable  to  attend. 

The  distinction  between  long  range  "normal  mode"  and  "other 
approaches"  is  somewhat  arbitrary  and  this  characterization  of  a  model 
was  made  by  AESD  prior  to  the  workshop.  The  NUSC-1  and  NUSC-2  models 
are  approximate  and  extend  normal -mode  models  respectively.  Figure  8 
lists  all  attendees  and  their  organizational  affiliations. 
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The  results  submitted  for  each  model  against  each  test  case  are 
presented  in  Figures  9  through  64.  As  a  reference  (but  by  no  means  a 
standard)  predictions  using  the  new  Navy  Interim  Standard  Transmission 
Loss  Model  (FACT)  have  been  included.  FACT  is  basically  a  ray  model 
which  includes  higher  order  wave  corrections  for  diffraction  and  sur¬ 
face-image  interference.  It  does  not  include  the  full  interference  of 
the  wave  field;  hence  predictions  of  transmission  loss  are  usually 
smoother  than  normal -mode  results. 


Predictions  For  Test  Case  1A  (Figures  9-20) 

The  predictions  for  Test  Case  1A  are  probably  the  most  consistent 
set  of  all  those  submitted.  Two  predictions  from  AESD  have  been  pro¬ 
vided.  The  first,  "ASRAP",  corresponds  to  the  surface-duct  module  in 
FACT  developed  by  Clay  and  Tatro  and  originally  implemented  in  RP70. 

It  is  basically  an  energy-conserving  model  intended  to  yield  averaged 
or  smoothed  ducted  levels.  The  second  prediction,  "FACT",  was  generated 
by  inhibiting  the  standard  Clay/Tatro  model,  tracing  rays  in  the  surface 
duct,  and  performing  the  standard  FACT  processing  of  these  rays.  When 
compared  with  the  normal -mode  results  quite  good  agreement  is  obtained, 
and  the  modal  interference  pattern  (only  2  modes  are  present)  is  seen 
to  correspond  to  the  small  convergence  zones  of  rays  in  the  duct.  It 
should  be  noted,  however,  that  for  the  case  selected  the  two  modes 
are  well  trapped  with  very  little  leakage  to  a  range  of  16nm.  Had  there 
been  more  leakage,  FACT  (which  cannot  include  leakage)  would  have  pre¬ 
dicted  too  high  an  intensity. 

Two  predictions  using  the  NUSC-3  model  were  provided,  one  allowing 
for  only  ducted  paths  and  another  allowing  for  the  bottom-reflected 
paths  present  for  the  full  sound-speed  profile.  The  second  has  been 
included  in  these  predictions  since  it  illustrates  the  deterioration  in 
the  smooth  modal  interference  pattern  expected  when  bottom-reflected 
energy  is  also  present. 


Predictions  For  Test  Case  IB  (Figures  21-32) 

The  primary  differences  between  the  various  predictions  for  this 
case  resulted  from  different  treatments  of  the  ocean  bottom.  This 
question  was  clearly  resolved  by  Bartberger  of  NADC  who  provided  three 
sets  of  predictions: 

1.  Modes  corresponding  to  refracted  and  RSR  paths  only  (with 
very  deep  nulls  between  the  convergence  zones). 

2.  These  modes  plus  those  corresponding  to  rays  reflecting 

at  less  than  the  critical  angle  (quite  similar  to  NUC-1 ,  BTL-1 ,  and  NOL). 

3.  These  modes  plus  the  contributions  from  modes  corresponding 
to  rays  striking  the  bottom  at  greater  then  critical  (NADC,  NUSC-1 , 
NUSC-2,  NUSC-3,  etc.)  In  the  convergence  zones  all  models  are  in  good 
agreement  except  perhaps  the  ARL  model  which  may  have  undersampled  the 
transmission-loss  curve. 
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Predictions  For  Test  Case  2  (Figures  33-43) 


The  same  level  of  agreement  is  present  for  this  case  as  for  case 
IB  with  differences  due  essentially  to  the  same  bottom-related  ques¬ 
tions.  The  agreement  between  FACT  and  the  mode  models  is  somewhat 
surprising  since  it  might  be  expected  that  transmission  to  a  receiver 
so  near  the  bottom  of  the  upper  duct  for  this  type  of  sound-speed  pro¬ 
file  would  be  inadequately  modeled  by  ray  theory. 


Predictions  For  Test  Case  3A  (Figures  44-54) 

In  this  case  several  modes  are  present  and,  with  the  exception  of 
the  ARL  model  (which  may  have  a  normalization-  error) ,  all  predictions 
are  in  good  agreement.  The  BTL-1  model  tends  to  show  slightly  higher 
losses  at  longer  ranges  than  the  mode  models,  and  the  source  of  this 
discrepency  is  not  resolved. 


Predictions  For  Test  Case 3B  (Figures  55-64) 

The  parameters  in  this  test  case  were  intentionally  selected  to 
correspond  to  a  duct  at  a  frequency  just  below  cut-off.  Unfortunately, 
the  hand  calculation  which  led  to  the  parameters  used  a  WKB  approxima¬ 
tion  to  determine  the  cut-off  and  the  consensus  was  that  probably  one, 
very  leaky  mode  was  present.  The  mode  predictions  reflect  this 
critical  regime  and  may  be  divided  into  approximately  three  categories: 

1.  Very  low  loss  at  range  (NUSC-2,  NUSC-4) 

2.  Moderate  loss  at  range  (NADC,  NRL-1 ,  NOL,  NRL-2) 

3.  High  loss  at  range  (ARL,  NUSC-1 ,  BTL-1). 

The  FACT  predictions  may  shed  some  light  on  these  differences. 

At  short  range  the  FACT  predictions  oscillate  rapidly  between  high- 
loss  bottom-bounce  regions  and  low-loss  RBR  regions.  Since  total 
phases  are  not  computed,  the  RBR  rays  propagate  successfully  and 
eventually  their  convergence  zones  overlap  so  that  the  FACT  prediction 
resembles  a  category  1  prediction.  At  short  ranges  where  the  zones 
are  very  narrow  FACT  resembles  the  type  2  category. 

The  large  differences  between  models  for  this  test  case  are  an 
indication  of  one  of  the  unresolved  areas  in  modeling.  It  was  also 
noted  that  the  introduction  of  a  very  slight  amount  of  absorption  in 
the  bottom  would  drive  the  category  2  predictions  down  to  those  in 
category  3.  The  sensitivity  of  the  predictions  to  such  parameters 
raises  several  additional  questions  for  those  involved  in  predictive 
acoustics. 
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BACKGROUND  FOR  ACOUSTIC  WAVE  MODELS* 


It  is  the  purpose  of  this  section  to  provide  the  common  points 
of  departure  for  the  various  classes  of  models  presented  at  the  work¬ 
shop.  The  emphasis  will  be  on  the  differences  between  these  classes 
with  the  subsequent  section  devoted  to  the  differences  between  specific 
models  within  a  class. 

For  a  medium  with  index  of  refraction 

n(r)  =  c0/c(r)  (1) 

where  c  is  sound  speed  and  c0  a  reference  sound  speed,  the  excess 
pressure,  P,  excited  by  a  harmonic  point  source  at  location  jr'  with 
angular  frequency  to  satisfies  the  inhomogeneous  reduced  wave  equation 

[V2  +  k2n2(rj]  P(r_,r/ )  =  -  6 (r-r.1 )  (2) 

where  the  wave  number  k  =  w/c0  and  r_  is  the  observation  point.  P  must 
satisfy  the  radiation  condition  at  infinity  and  vanish  at  the  ocean  surface. 

At  this  point  nearly  all  models  assume  that  the  ocean  depth  and 
sound  velocity  profile  are  independent  of  range.  Then  (2)  can  be 
expressed  in  cylindrical  coordination  as 

p  &P&  +  1F+  k2,,2(z)]  p(^'>  '  "Zip  «(»>«<*-*'>•  <3> 

The  solution  of  this  boundary  value  problem  can  be  represented  in 
many  forms,  each  having  different  properties.  Four  representations  are 
described  here  which  form  the  basis  for  most  of  the  models  presented. 

The  first  is  a  contour-integral  representation  in  the  plane  of  the 
separation  constant.  Labianca  has  shown  that  the  other  three  (Fourier- 
Bessel  integral,  the  residue  series,  and  the  normal-mode  solution)  are 
derivable  from  the  X-plane  representation. 


The  X-Plane  Representation 

The  X-plane  representation  is  a  contour-integral  representation 
which  is  applicable  to  any  separable,  cylindrical -coordinate,  boundary- 
value  problem  with  symmetry  about  the  z-axis.  For  the  case  of  equation 
(3),  the  solution  takes  the  following  form: 


*  The  following  discussion  has  been  extracted  from  "Normal  modes, 
virtual  modes  and  alternative  representations  in  the  theory  <\f  surface- 
ducted  sound  propagation"  by  F.M.  Labianca,  JASA  53  (1973),  1137-1147. 
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G-j  (p;l-X)G2(z,z'  ;X)dX  . 


(4) 


P(r,r') 


Here,  G^(p;l-X)  and  G2(z,z';X)  are  the  one-dimensional  Green's  functions 

for  the  p  and  z  coordinates,  respectively.  They  satisfy  the  differential 
equations 


1  d 

P  3p  p 


d 


+  r(i-x) 


G-j  (p;l-X) 


(5a) 


and  (letting  q(z)  =  l-n2(z)) 
„  1 


+  k2(X-q(z) ) 
dz 


G2(z,z';X) 


«(z-z') 


(5b) 


defined,  respectively,  on  the  domains  0  £  p  <_  p^  and  ±z  <_  z 2*  It  is 

seen  from  Eq.  5a  that  G^ (p ;  1  -X)  is  required  to  satisfy  the  delta-function 

source  condition  at  p  =  0,  which  is  one  of  the  endpoints  of  its  domain. 

In  addition,  it  is  required  to  satisfy  a  linear,  homogeneous  boundary 
condition  at  the  other  endpoint  p^ .  On  the  other  hand,  G2(z,z';X) 

satisfies  the  source  condition  at  a  point  z'  within  the  endpoints  of  its 
domain  and  is  required  to  satisfy  linear,  homogeneous  boundary  conditions 
at  the  endpoints  z^  and  z2> 


The  contour  C  ,  over  which  the  integral  in  Eq.  4  is  to  be  evaluated, 

X 

encloses  the  X-plane  singularities  of  G2  in  a  positive  sense  (counter¬ 
clockwise)  or  the  singularities  of  G^  in  a  negative  sense  (clockwise). 

G^  is  given  by 

G-j  (p ;  1  -X)  =  1  h£1}  (k p),  (6) 

being  the  Hankel  function  of  the  zeroth  order  and  first  kind.  In 

order  to  insure  that  the  radiation  condition  is  satisfied  for  all  complex 
values  of  X,  the  function  /I -X  is  made  single-valued  by  cutting  the 
X-plane  from  -«  to  1  and  requiring  /1-X  >  0  on  the  entire  top  Riemann 
sheet.  The  contour  is  then  taken  to  encircle  this  cut  in  a  negative 
sense. 

The  function  G2(z,z';X)  is  the  characteristic  Green's  function  of 
the  depth  geometry,  and  it  is  required  to  vanish  at  the  endpoint  z^  =  0 
and  to  satisfy  the  radiation  condition  at  the  endpoint  z^  =  Note  also 
that  G2(z,z';X)  and  dG2(z,z' ;X)/dz  must  be  continuous  for  all  z.  The 
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solutions  for  the  eigenvalues,  X,  and  G2  basically  differentiate  the 
types  of  models  within  each  class. 

The  Fourier-Bessel  Representation 

The  Fourier-Bessel  representation  is  obtained  from  the  X-plane  repre¬ 
sentation  by  a  change  of  variables  and  a  contour  deformation  to  yield 

00 

p(r,r')  =  2~  l  cJ0(cp)G2(z,z' ;l-c2/k2)dc  .  (7) 

J  0 

The  Residue-Series  Representation 

The  contribution  to  the  integral  from  the  poles  within  a  deformed 
contour  gives  the  solution 

oo 

P(r.r')  =  ^  JfH0^(k  /l_Aj  P^jU(z,X_.  )U(z' ,X. )  (8) 

J  1 

where  the  U's  play  the  role  of  eigenfunctions  and  the  Xj  must  be  found 

in  the  complex  X  plane.  The  contour  deformation  which  yields  this  re¬ 
sult  is  val id  only  for 

P  >  J?  (z+z1 )  ;  (9) 

hence  the  Uj's  are  not  eigenfunctions  (in  the  sense  that  they  do  not 
form  a  complete  set)  and  the  summation  diverges  for  p  <  ^T(z+z'), 


The  Normal  Mode  Representation 

A  different  contour  deformation  leads  to  the  normal  mode  representation 


r°° 

P(r,r' )  =  \  (k  /NX  p)|i(z,X)i(z' ,X)dX, 

J  -00 


(10) 


where  $  are  related  to  the  true  eigenfunctions  <J>. 

If  at  some  depth  the  ocean  bottom  Is  considered  to  be  Isovelocity 
with  a  sound  speed  greater  than  the  minimum  sound  speed  in  the  water  then 
some  of  the  Xj  are  real  valued  and  the  integral  in  Eq.  10  becomes  a  sum¬ 
mation  over  these  discrete  modes  plus  an  integral  over  the  continuous 
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spectrum  corresponding  to  the  remaining  set  of  complex  X.,  In  terms  of 

J 

a  ray  description,  the  discrete  modes  correspond  to  rays  which  are 
either  refracted  in  the  water  column  or  strike  the  bottom  at  less  than 
the  critical  angle.  The  continuous  spectrum  (and  complex  X-)  corresponds 

\J 

to  partially  reflected  rays  striking  the  bottom  at  greater  than  critical 
plus  diffraction  fields  such  as  "head  waves". 


Classification  of  Models 


With  this  background  and  from  the  presentations  of  the  participants 
we  have  attempted  to  classify  each  of  the  wave  models: 

1.  X-plane  representation  -  ARL 

2.  Fourier-Bessel  representation  -  NUSC-3 

3.  Residue  series  -  NUC-1 ,  ARL,  NADC  (special  purpose  models), 
CU,  NRL-2,  NUC-2,  LDGO 

4.  Normal  mode  -  NRL-1 ,  NADC  (general  ocean  model),  APL,  NOL, 
NUSC-4,  NUSC-1 ,  NUSC-2 ,  NRL-3,  BTL-2 

As  can  be  seen  the  initial  classification  of  models  in  Figure  7  did 
not  distinguish  between  normal  modes  and  residues, and  the  NUSC-1  and 
NUSC-2  models  were  put  in  the  "Other  Approaches"  category  since  they 
were  rather  different  from  the  other  mode  programs. 

The  BTL-1  model  of  Tappert  and  Hardin  does  not  fall  into  any  of  the 
above  categories  since  it  departs  from  all  other  models  by  approximating 
the  wave  equation  and  numerically  solving  the  resulting  equation  without 
resort  to  separation  of  variables  or  normal  modes.  This  approach  will 
be  discussed  in  detail  in  the  section  devoted  to  it. 


PRESENTATION  OF  MODELS 


This  section  summarizes  the  presentations  of  each  of  the  models 
emphasizing  the  differences  between  them.  Any  inaccuracies  in  these 
descriptions  are  regretted  and  readers  are  advised  to  contact  the 
principal  investigators  for  detailed  descriptions  of  their  models. 


NUC-1  -  Naval  Undersea  Center  -  D.  Gordon  and  M.  Pedersen 

This  is  a  residue  model  in  which  the  sound  speed  profile  (range- 
independent)  is  segmented  into  an  arbitrary  number  of  layers  within  which 
n2(z)  varies  linearly.  This  leads  to  a  modified  equation  for  the 
residues  X,  and  pseudo-eigenfunctions  U  of  the  form 

gp-  +  ^k2-X2  -  a^z-z^U  =  0  (11) 
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where  z^  is  the  layer  depth  and  a^  depends  on  the  sound-speed  and 

dn2/dz  in  the  layer.  The  technique  assumes  an  implicit  transformation 
such  that  this  depth-dependent  equation  is  a  standard  form.  Either, 
for  a  given  c(z),  the  independent  variables  are  transformed  such  that 
the  resultant  equation  is  of  a  standard  form  or,  given  the  equation 
with  the  independent  variable  being  z,  the  speed  of  sound  c(z)  is 
transformed.  The  deepest  layer  is  a  semi-infinite  space  with  a  nega¬ 
tive  sound-speed  gradient.  At  the  bottom,  a  radiation  condition  is 
imposed  such  that  only  outgoing  waves  are  retained. 

In  order  to  determine  which  mode  is  in  fact  being  calculated'  a 
transformation  is  employed  which  essentially  corresponds  to  the  Airy 
function  behavior  of  the  depth  functions  in  each  layer.  This  allows 
a  check  not  only  on  which  mode  is  being  calculated  but  also  that  all 
modes  have  been  found.  The  model  is  not  particularly  sensitive  to 
the  choice  of  layer  depths  or  the  number  of  layers.  A  complex  sound 
speed  may  be  specified  in  each  layer  to  simulate  volume  absorption, 
though  generally  this  is  done  only  for  the  bottom. 

The  time  and  cost  required  for  the  test  cases  ib  a  UNIVAC  1108  were 
as  follows: 

Cases  1A,  IB:  36  sec.  CPU  ($6) 

Case  2:  (120  modes)  1  min.  56  sec.  ($19) 

Case  2  Contours  of  loss  in  range  and  depth:  3  min,  43  sec.  ($36). 

The  program  is  written  in  double  precision  -  with  4,000  program  and 
1,700  library-subroutine  steps.  There  is  a  50K  core  storage  limit  which 
implies  350  modes,  8  layers  and  50  receiver  depths. 


NRL-1  -  Naval  Research  Laboratory  -  J.  C.ybulski 

This  is  the  oldest  model  and  is  the  basis  for  a  number  of  normal¬ 
mode  models  presented  by  the  other  Navy  laboratories.  The  model  assumes 
a  range-independent  environment  in  which  the  sound  speed  varies  arbitrarily 
with  depth  and  is  constant  in  the  bottom.  Finite-difference  techniques 
are  employed  using  a  numerical  integration  of  the  depth-dependent  solutions. 
An  initial  guess  for  the  eigenvalue  is  obtained  by  assuming  the  speed  of 
sound  in  water  is  constant.  The  integration  proceeds  from  the  ocean 
bottom  to  the  surface  and  if  the  boundary  condition  P(0)  =  0  is  not  sat¬ 
isfied  the  eigenvalue  is  altered  slightly  and  the  process  repeated.  For 
each  choice  of  eigenvalue,  the  zero  crossings  are  counted  assuring  that 
no  modes  are  missed.  As  with  most  models  the  solution  employs  the 
asymptotic  form  of  the  Hankel  function  valid  at  distances  many  wave¬ 
lengths  from  the  source. 
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For  Case  IB,  43  modes  were  required  taking  6  minutes  on  a  CDC  3800; 
for  Case  2,  50  modes;  for  Case  3A,  6  modes;  and  for  Case  3B  one  mode. 
The  core  storage  for  the  program  is  73K  octal.  All  the  eigenvalues  are 
real.  (If  complex  eigenvalues  are  generated  they  are  discarded.)  For 
example,  only  one  real  eigenvalue  was  found  for  Case  3B,  corresponding 
to  one  trapped  mode. 


NADC  -  Naval  Air  Development  Center  -  C.  Bartberger 

Several  NADC  models  have  been  developed  and  can  be  applied 
under  a  variety  of  conditions.  Restricted  models  treat  "oceans"  with 
no  bottom,  and  1,  2  or  3  layers  which  are  special  cases  of  the  NUC 
Residue  model.  The  square  of  the  index  of  refraction  varies  linearly 
in  each  layer  yielding  Airy  functions  for  the  eigenfunctions.  The 
bottom,  when  multi-layered,  allows  for  either  volume  attenuation  or  a 
bottom  loss  specified  as  a  function  of  the  grazing  angle  for  the  ray- 
equivalent  of  the  mode. 

A  complete  ocean  model  has  also  been  developed  which  allows  for  50 
layers  in  the  sound-speed  profile.  The  original  input  profile  is  then 
effectively  approximated  by  homogeneous  layers  in  a  finite-difference 
algorithm  similar  to  that  in  NRL-1 .  Initial  guesses  for  the  eigenvalues 
are  obtained  using  the  WKB  approximation  and  the  solution  is  obtainedbfcy 
integrating  in  depth  from  the  bottom  to  the  surface,  counting  zero 
crossings  and  iterating.  While  the  models  are  fully  automated  they 
sometimes  have  difficulty  finding  the  eigenvalues.  Bartberger  feels  that 
the  arbitrary  profile  model  is  as  flexible  as  the  present  ray  models. 

He  suggests  that  an  obvious  improvement  is  to  have  a  variable  layer 
thickness.  Furthermore,  he  believes  that  there  should  be  at  least  two 
models:  a  surface-duct,  bottomless  model  and  a  full -ocean  model  incor¬ 
porating  bottom  characteristics. 


ARL  -  Applied  Research  Laboratory  -  R.  Deavenport  and  J.  Beard 

The  ARL  model  uses  the  A-plane  representation  for  a  specialized 
case  where  the  inverse  square  of  the  velocity  field  is  characterized  by 
a  five  parameter  function,  the  Epstein  profile: 

p  =  A  sech2  ^p-  +  B  tanh  ^jp-  +  D  (12) 

This  is  extendable  to  seven  parameters  and  the  coefficients  may  be  deter¬ 
mined  to  reasonably  model  refracted  energy  in  deep-ocean  sound  channels. 
Surface  or  bottom-reflected  paths  cannot  be  modeled.  The  depth-dependent 
Green's  function  (z ,z 1 ;d )  is  given  by  hypergeometric  functions  and  the 

eigenvalue  spectrum  is  composed  of  discrete  values  and  a  continuous  set. 
No  approximations  are  employed  in  arriving  at  the  solution,  however,  it 
has  been  found  to  be  very  sensitive  to  small  changes  in  the  velocity 
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profile.  The  program  is  very  rapid  and  was  presented  as  a  possible 
control  program  to  test  approximate  techniques.  The  solutions,  however, 
only  include  those  modes  with  phase  velocities  close  to  the  sound  speed 
at  the  source. 

There  was  substantial  discussion  regarding  the  ARL  output  and  its 
relation  to  the  NUC  output  at  short  ranges  for  Case  IB  since  the  "bottom¬ 
less"  ARL  model  showed  higher  levels  between  10  and  20  miles  than  the 
NUC  model  with  bottom-reflected  paths.  In  the  course  of  the  discussion 
it  was  noted  that,  in  bottom-bounce  areas  with  high  loss,  important 
contributions  to  the  field  come  from  only  a  few  modes  which  are  always 
in  phase;  the  remaining  modes  tend  to  destructively  interfere  with  each 
other  and  the  resulting  low  level  may  have  been  missed  in  the  ARL  model 
by  limiting  the  modes  considered. 

APL  -  Applied  Physics  Laboratory  -  N.  Nicholas 

The  APL  model  is  also  a  sound-channel  model,  (no  surface  or 
bottom)  in  which  the  sound-speed  profile  is  parabolic  in  depth: 

cmm  0  '  “2z2),/2  031 

The  depth  functions  for  the  resulting  normal  modes  are  Hermite  Poly¬ 
nomials  and  are  finite  in  number  for  this  type  of  profile.  This  model 
has  been  used  primarily  to  investigate  the  periodic  focusing  regions 
predicted  for  parabolic  profiles  in  axis-to-axis  geometries.  Slightly 
differing  velocity  profiles  may  be  accomodated  using  perturbation 
techniques  analogous  to  those  of  quantum  mechanics.  The  corresponding 
potential  function  is  of  the  form 


V  (z)  =  k2ct2z2  +  V:(  z)  (14) 

where 

\l1(z)  «  k2a2z2  (15) 

and  V*(z)  is  a  fourth  order  polynomial.  Strong  axial  focusing  disappears 
with  the  introduction  of  the  perturbation  potential.  Modifications  to 
include  boundaries  have  been  explored  and  give  rise  to  infinite  series 
in  addition  to  the  Hermite  Polynomials.  The  radius  of  convergence  of 
these  series  of  solutions  is  quite  reasonable  when  tested  with  conven¬ 
tional  power-series  convergence  criteria. 


NOL  -  Naval  Ordnance  Laboratory  -  I.  Blatstein,  H.  Uberall ,  A.  Newman 

The  NOL  model  is  an  outgrowth  of  the  NRL-1  normal -mode  model 
applicable  to  deep  water,  and  uses  a  finite-difference  technique  to 


obtain  numerical  solutions.  The  continuous  portion  of  the  spectrum  is 
ignored  for  deep  water  cases  and  only  modes  with  phase  velocities  less 
than  the  velocity  of  the  bottom  are  included.  The  equivalence  of  the 
normal -mode  and  residue  techniques  is  demonstrated  by  the  excellent 
agreement  with  NUC-1  for  Case  IB.  The  following  table  summarizes 
running  statistics  of  the  NOL  model  for  the  four  test  cases: 


Case 

Modes 

Time  (sec) 

Cost 

IB 

43 

230 

$(25.00)? 

2 

81 

350 

$  40.00 

3A 

5 

19 

$  2.40 

3B 

1 

42 

$  4.70 

All  numerical  solutions  are  single  precision  and  core  storage  is  60K 
octal  words  on  a  CDC  6400. 

Comparisons  were  shown  of  predicted  convergence-zone  levels  using 
ray  theory,  normal-mode  and  modified  ray  theory  as  given  by  Silbiger 
and  Sachs.*  Within  the  insonified  region  the  three  approaches  are  in 
very  good  agreement,  until,  in  the  immediate  vicinity  of  the  caustic, 
the  ray  theory  intensity  becomes  infinite.  In  the  shadow  zone  the 
normal-mode  and  modified  ray-theory  results  maintain  a  high  level  of 
agreement. 


CU  -  Catholic  University  -  V.  Nomady  and  H.  liberal! 

This  program  was  specifically  designed  to  consider  the  effects 
of  earth  curvature  on  wave  solutions  for  a  deep  ocean  at  long  ranges. 
The  ocean  is  characterized  by  a  single  layer  in  which  the  sound  speed 
is  given  by  an  Epstein  profile.  The  field  is  then  expressed  in  terms 
of  spherical  harmonics  (Legendre  and  Jacobi  polynomiols).  By  applying 
the  Watson  transformation  the  CW  response  is  obtained  in  terms  of  two 
contour  integrals  plus  a  series  of  residues  where  the  series  dominates 
at  long  ranges.  Physically  the  solution  may  be  interpreted  as  standing 
waves  in  the  depth  coordinate  and  circumferentially  outgoing  waves  on  a 
spherical  boundary.  Efforts  are  under  way  to  obtain  the  solution  for  a 
pulsed  source. 


NRL-2  -  Naval  Research  Laboratory  -  R.  Fitzgerald 

This  model  is  essentially  the  same  as  the  NUC-1  model;  the  major 
differences  being  associated  with  the  treatment  of  the  ocean  bottom. 

The  NUC-1  model  allows  the  sound  speed  to  go  to  zero  for  very  large  z, 
whereas  this  model  terminates  the  profile  exactly  as  specified  allowing 
for  fully  trapped  modes.  As  presently  configured  the  program  does  not 
attempt  to  normalize  the  source.  The  results  were  normalized  by  AESD 
by  comparing  key  features  in  the  predictions  with  similar  features  in 
other,  normalized,  predictions. 


*  JASA,  Volume  49,  1971,  pp.  824-840 
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NUSC-4  -  Naval  Undersea  Systems  Center  -  W.  Kanabis 


This  program  Is  a  normal-mode  model  nearly  identical  to  NRL-1 
and  was  designed  primarily  for  shallow-water  studies.  The  only  modes 
which  are  kept  correspond  to  bottom  grazing  angles  which  are  less  than 
critical  and  these  modes  suffer  no  loss  due  to  interaction  with  the 
bottom.  For  current  applications  no  normalization  is  attempted. 


NUSC-1  -  Naval  Undersea  Systems  Center  -  G.  Leibiger 

This  model  is  called  "Ray -Mode"  since  it  uses  ray  concepts  to 
assist  in  obtaining  a  wave  solution.  Expressing  the  field  in  a  Green's 
function  formulation  the  zeros  of  the  Wronskian  (or  poles  of  the  Green's 
function)  are  found  by  phase-integral  techniques  using  rays.  From 
this  WKB  technique  the  real  and  imaginary  parts  of  the  eigenvalues  are 
determined.  Since  the  X-plane  integral  is  dominated  by  contributions 
near  these  poles  the  integral  may  be  replaced  by  a  summation  and  the 
contributions  from  each  term  may  be  interpreted  physically  in  terms  of 
rays  cycling  through  the  medium.  While  the  model  was  presented  as  it 
might  apply  to  a  surface-duct  problem,  the  technique  is  generally 
applicable  by  an  appropriate  modification  to  the  Wronskian. 


NUSC-2  -  Naval  Undersea  Systems  Center  -  H.  Weinberg 

This  is  a  normal -mode  model  capable  of  treating  a  fully  three- 
dimensional  environment  (that  is  both  range  and  azimuth  dependences  in 
the  sound-speed  and  bottom  descriptions).  By  employing  a  stretching  of 
coordinates  in  the  horizontal  plane  with  a  small -parameter  expansion 
the  field  is  described  in  terms  of  the  normal -modes  or  depth  functions 
of  the  geometry  which  areaallowed  to  vary  slowly  with  range.  (In  the 
present  version  only  adiabatic  variations  are  permitted.)  Departures 
from  cylindrical  symmetry  are  accounted  for  by  computations  of  the 
azimuthal  divergence  in  the  amplitude  of  each  mode  as  given  by  a  ray- 
trace  in  a  coordinate  system  orthogonal  to  the  modes.  In  practice  the 
modes  must  be  recomputed  at  each  range  where  the  sound-speed  profile 
or  the  bottom  depth  changes  significantly.  Predictions  were  compared 
with  measurements  made  in  a  range-dependent  environment  showing  excellent 
agreement  to  very  long  ranges. 


NUSC-3  -  Naval  Undersea  Systems  Center  -  F.  DiNapoli 

This  fast  field  program  (FFP)  directly  integrates  the  Fourier- 
Bessel  transform  allowing  for  a  segmented  (range  independent)  velocity 
profile  and  a  fluid  bottom.  After  taking  the  asymptotic  expansion  for 
Jg,  the  integral  was  found  to  be  strikingly  similar  to  a  Fourier  Trans¬ 
form.  Thus,  the  solution  appeared  to  be  a  discrete  Fourier  Transform 
from  wave  number  space  to  range  space  where  the  Green's  function  was 
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considered  to  now  be  a  function  of  wave  number.  The  difficulty  was 
to  obtain  a  Green's  function  which  permitted  rapid  computation  for 
many  wave  numbers.  Various  techniques  have  been  tried  and  abandoned 
including  numerical  integration  and  solutions  in  terms  of  hypergeometric 
functions.  Known  solutions  were  then  investigated  and  it  was  found  that 
if  the  sound  velocity  were  described  in  terms  of  exponentials  then  the 
Green's  function  became  a  linear  combination  of  Beisel  functions.  In 
fact,  the  wave  number  appears  only  in  the  order  of  the  Bessel  function 
and  recursion  relations  were  found  which  made  the  computation  feasible. 
In  practice  the  input  sound-velocity  profile  is  fit  with  several  ex¬ 
ponential  segments  which  then  form  the  working  profile.  The  program  is 
three  years  old  and  has  yielded  excellent  agreement  with  NADC's  model 
on  several  test  cases.  Running  times  are  comparable  to  those  for  the 
NADC  and  NOL  models  and  the  frequency  may  range  from  very  low  to  perhaps 
3.5kHz.  This  is  the  only  model  which  can  routinely  accomodate  such  high 
frequencies.  Such  cases  are,  however,  severely  limited  in  range. 


LDGO  -  Lamont-Dougherty  Geophysical  Observatory  -  H.  Kutschale 

This  program  is  primarily  concerned  with  under-ice  acoustics  and, 
therefore,  allows  for  a  solid  surface  as  well  as  a  solid  bottom.  The 
solution  consists  of  a  branch-line  integral  plus  residue  contributions, 
and  FFP  techniques  (see  NUSC  Rept.  No.  4103)  have  been  incorporated  to 
evaluate  the  branch-line  integral.  A  liquid  bottom  gives  rise  to  one 
branch-line  integral  while  a  solid  bottom  gives  rise  to  two.  The  con¬ 
tributions  to  the  field  from  the  branch-line  integral  are  felt  to  be 
quite  important  at  low  frequencies. 


BTL-1  -  Bell  Telephone  Laboratories  -  F.  Tappert  and  R.  Hardin 

This  model's  last  point  of  commonality  with  all  other  models 
presented  is  the  wave  equation.  The  model  has  three  basic  components, 
each  of  them  independent  of  the  other  two: 

1.  The  model  approximates  the  elliptic  wave  equation  with 
a  parabolic  equation  equivalent  to  a  small-angle  approximation  in  ray 
theory. 

2.  The  resulting  equation  is  then  solved  by  means  of  the 
recently  developed  split-step  Fourier  algorithm. 

3.  The  resulting  levels  of  intensity  as  a  function  of  range 
and  depth  are  then  contoured  using  a  new,  rapid  contouring  technique. 

Beginning  with  the  wave  equation  for  a  cyl indr i cal ly  symmetric  geo¬ 
metry  and  letting 

P(r,z)  =  ip(r ,z)  hJ2^  (kr)  (16) 

and  making  the  following  approximations: 
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kr  »  1 


Vrr  «  2i  kipr  (17) 

the  parabolic  equation 

i>r  +  2F^zz  +  |  (n1 2 3 4(r,z)-l )  iff  =  0  (18) 

is  obtained.  The  application  of  this  parabolic  equation  to  a  number  of 
problems  is  discussed  in  Electromagnetic  Diffraction  and  Propagation 
Problems  edited  by  Fock  (Peagamon  Press,  1965).  The  approximations  which 
yield  the  parabolic  equation  are:  first,  the  usual  far-field  asymptotic 
expansion  for  the  Hankel  function;  and  second,  the  assumption  that  the 
waves  are  predominantly  radial  -  that  is  the  corresponding  rays  are  at 
relatively  shallow  angles.  This  second  approximation  includes  the 
neglect  of  any  backscatter. 

For  the  original  elliptic  wave  equation  the  boundary  conditions 
must  be  prescribed  on  a  closed  surface  (or  as  a  radiation  condition), 
and  the  direct  solution  involves  iterative  approaches  such  as  relaxation 
techniques.  Parabolic  equations  may  be  solved  by  marching  the  solution 
away  from  a  surface  on  which  initial  conditions  are  specified.  Hence, 
if  ip(r,z)  is  specified  for  any  range  r  then  the  split-step  Fourier 
algorithm  yields  at  a  small  step  Ar: 


i^(r+Ar,z)= 


e14r  k  <"’-!>/*  3- V^/a9Wr(2))) 


(19) 


where  ,,is  a  Fourier  transform  from  z  to  i  and  is  the  inverse 

transform  from  £  to  z. 

Application  of  this  technique  to  ocean-acoustics  problems  requires 
the  following  steps: 

1.  Define  a  basic  medium  which  consists  of  the  original  medium 
reflected  across  the  ocean  surface. 

2.  Insert  into  this  geometry  the  source  (or  receiver)  of  interest 
and  its  reflection  (or  image)  phase  shifted  by  180  degrees.  This  removes 
the  boundary  condition  at  the  ocean  surface. 

3.  Extend  the  medium  (and  its  reflection)  with  a  highly  absorbing 
layer  at  a  depth  greater  than  the  deepest  portion  of  the  ocean  considered. 

4.  This  double  medium  may  now  be  considered  to  be  infinitely 
periodic  in  depth  since  no  energy  from  the  other  pairs  of  sources  can 
penetrate  the  absorbing  layer  to  affect  the  field  at  points  of  interest. 
Hence,  the  region  of  interest  may  be  sampled  with  a  grid  in  z  and  Fast 
Fourier  Transforms  may  be  used  to  evaluate  Eq.  19  resulting  in*a  rapid 
integration  of  \p. 


15 


The  graphical  contour  displays  are  quite  useful  in  analyzing  the 
field  properties  and  examples  were  shown  for  each  of  the  test  cases 
plus  two  range-dependent  problems: 

1.  Long-range  transmission  in  the  North  Pacific. 

2.  The  field  in  a  weakening  surface  duct.  Figures  65  and  66 
illustrate  this  case  for  the  surface  duct  of  test  case  1A  which  is  then 
allowed  to  become  isovelocity  (the  range  scale  of  this  second  plot  is 
slightly  more  compressed  than  on  the  previous  plot). 

Tappert  outlined  the  strengths  and  weaknesses  of  the  model.  The 
obvious  strengths  in  terms  of  the  environment  are: 

1.  A  range-dependent  sound-speed  profile. 

2.  A  variable-depth  ocean  bottom  (which  is  modeled  as  refracting 
rather  than  as  a  reflecting  surface). 

3.  Variable  volume  absorption  may  be  included  directly  in  the 
index  of  refraction. 

The  limitations  of  the  model  are: 

1.  The  small  angle  approximation  limiting  equivalent  rays  to 

=  20°. 

2.  The  resulting  limitation  that  energy  reflecting  from  steep 
up-slopes  into  steeper  angles  may  be  improperly  modeled. 

3.  The  high-frequency  limit  determined  by  the  maximum  size  FFT 
which  can  efficiently  be  used. 

4.  Discontinuities  requiring  careful  treatment  to  avoid  numerical 
problems. 

5.  The  present  model  was  not  normalized  and  there  are  some  prob¬ 
lems  associated  with  modeling  a  point  source.  AESD  normalized  the  pre¬ 
dictions  in  the  same  way  as  the  NRL-2  predictions. 

The  advantages  of  the  numerical  technique  which  arise  in  the  solution 
of  parabolic  equations  are: 

1.  The  accuracy  is  good  to  second  order  in  range  which  implies 
large  steps  in  range  are  possible. 

2.  Energy  is  conserved. 

3.  The  numerical  technique  is  unconditionally  stable. 

4.  It  is  easy  to  implement. 

5.  It  is  computationally  efficient. 

The  disadvantages  are: 

1.  The  boundary  conditions  must  be  periodic  in  depth. 

2.  A  uniform  grid  is  required. 

A  number  of  questions  remain  concerning  specific  applications  to  ocean 
acoustics  problems  since  the  program  was  originally  developed  for  electro¬ 
magnetic  applications.  The  consensus  at  the  close  of  the  workshop  was, 
however,  that  this  was  one  of  the  most  exciting  innovations  in  propagation 
modeling  in  recent  years. 
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NRL-3  -  Naval  Research  Laboratory  -  F,  Ingenito 


This  is  a  normal -mode  program  similar  to  the  NRL-1  and  NOL 
models  except  that  it  is  limited  to  three  layers.  It  calculates  only 
discrete  modes  and  has  real  eigenvalues  and  real  eigenfunctions.  At 
present  it  is  being  expanded  to  consider  bottom  and  surface  absorption. 


NUC-2  -  Naval  Undersea  Center  -  H.  Morris  and  H.  Bucker 


This  is  the  same  basic  residue  model  as  NUC-1  but  restricted 
to  surface-duct  geometries  and  extended  to  include  rough-surface 
effects.  Lambert's  Law  is  applied  to  the  modes  (with  the  angle  of 
incidence  obtained  from  the  ray  equivalent)  to  describe  rough-surface 
scattering.  The  first-order  effects  are  a  reduction  in  intensity  in 
the  duct  and  an  enhancement  below  the  duct.  Second  order  effects  are 
being  investigated  which  include  cross-coupling  between  modes.  Com¬ 
parisons  with  SUDS  data  showed  impressive  agreement  in  many  cases 
though  the  surface-duct  problem  is  by  no  means  completely  solved.  This 
model  is  described  in  detail  in  JASA,  volume  48,  1970,  pp.  1187-1194. 


BTL-2  -  Bell  Telephone  Laboratories  -  F.  Labianca  and  E.  Harper 

This  is  also  a  surface  duct  model  using  normal  modes  rather 
than  residues.  The  basic  model  is  documented  in  JASA  volume  53,  1973, 
pp.  1137-1147.  Because  normal  modes  are  used,  rather  than  residues,  a 
basis  is  available  for  expanding  arbitrary  functions.  This  allows  for 
a  novel  approach  to  the  incorporation  of  a  rough  surface  in  the  surface 
duct  by: 


1.  Decomposing  the  surface  into  its  Fourier  components 
(assuming  the  applicability  of  the  method  of  superposition,  i.e.  that 
the  surface  may  be  linearized). 

2.  Expanding  the  surface  profile  for  each  component  in 
terms  of  the  eigenfunctions  or  normal  modes  of  the  surface  duct. 

3.  Considering  realistic  surfaces  with  small  displacements 

and  slopes. 

4.  Using  perturbation  techniques  on  the  basic  mode  solution 
to  compute  the  resulting  rough-surface  effects. 

Questions  raised  and  not  fully  answered  about  this  technique  concern 
the  domains  of  applicability  in  terms  of  the  relative  size  of  acoustic 
wavelength,  surface  wavelength,  surface  amplitude,  and  duct  thickness. 


WORKING  SESSIONS 


It  was  the  purpose  of  the  working  sessions  to  focus  the  discussion 
on  specific  applications  and  to  assess  the  state-of-the-art  in  these 
areas.  While  issues  are  rarely  settled  at  such  sessions  (and  frequently 
more  heat  than  light  is  generated)  the  needs  of  the  advanced  development 
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community,  as  expressed  by  AESD,  were  conveyed  to  the  participants.  It 
is  impossible  to  present  fully  the  discussions  in  each  working  session 
and  in  this  report  we  will  only  attempt  to  reproduce  the  basic  questions 
and  our  resulting  assessment  of  the  state-of-the-art. 


A.  Rough-Surface  Scattering 

The  basic  question  is:  how  may  a  rough  ocean  surface  be  incorporated 
in  acoustic  wave  propagation  programs. 

The  NUC-2  and  NUSC-1  models  are  the  only  ones  which  currently  con¬ 
sider  rough-surface  effects.  For  the  coherent  field  an  effective 
reflection  coefficient  is  used.  The  incoherent  field  is  scattered  in 
all  directions  using  Lambert's  Law.  The  mode  reflection  coefficient  is 
inferred  from  the  plane-wave  reflection  coefficient  for  an  angle  of 
incidence  given  by  the  ray  equivalent  of  the  mode  and  is  limited  to 
statistically  rough  surfaces  of  small  amplitude.  The  parameter  required 
is  "significant  wave  height"  which  may  be  related  to  rms  wave  amplitude 
and  thence  to  wind  speed. 

The  BTL -rough-surface  approach  described  by  Labianca  is  the  only  one 
under  consideration  which  does  not  require  the  mixture  of  ray  and  mode 
concepts  and  hence  is  on  a  rigorous  theoretical  basis.  The  approach 
has  not  yet  been  implemented  and  the  questions  regarding  relative  lengths 
remain  to  be  resolved. 

A  detailed  presentation  of  the  mathematical  formulation  of  surface 
scattering  was  made  by  Richard  Hoi  ford  of  BTL.  He  presented  an  exact 
analytical  solution  for  plane-wave  scattering  from  a  periodic  surface. 

The  basic  assumption  of  physical  optics  is  not  made  and  the  solution  is 
presented  in  the  form  of  an  integral  equation  of  the  second  kind  which 
is  solved  by  matrix  inversion  where  convergence  is  assured.  This 
technique  may  serve  as  a  useful  control  for  comparisons  with  approximate 
techniques  and  Hoi  ford  showed  a  sample  case  where  the  method  of  physical 
optics  clearly  breaks  down.  This  approach  to  scattering  is  available  in 
several  BTL  reports. 

George  Zipfel  presented  work  from  a  paper  he  and  John  Di santo  of 
NRL  published  on  surface  scattering  (Journal  of  Math  Physics  13,  1972, 
pp.  1903.)  The  surface  is  treated  as  a  multi-variant  Gaussian  distribu¬ 
tion  and  the  results  are  presented  as  an  ensemble  average  in  terms  of 
correlation  functions.  A  series  of  diagramatical  rules  similar  to 
Feynman  diagrams  enable  a  straightforward  approach  to  the  problem  to  be 
implemented. 

With  the  above  techniques  the  modal  approaches  to  incorporation  of 
rough  surface  in  wave  programs  appeared  to  be  heading  in  the  right 
direction  for  the  problem  of  paramount  importance  -  the  surface  duct. 
Tappert  felt  that  the  parabolic  equation  technique  could  be  applied  to 
rough  surfaces,  especially  for  the  surface  duct  where  the  angles  of 
incidence  are  small. 
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B.  Volume  Attenuation 


The  volume  attenuation  question  is  primarily  one  of  measurement. 

While  measuring  volume  attenuation  at  high  frequencies  is  relatively 
straightforward,  low-frequency  measurements  require  extremely  long 
baselines  and  have  led  to  several  conflicting  results.  The  basic 
technique  is  to  measure  long-range  transmission  loss  at  a  number  of 
frequencies  in  an  axis-to-axis  geometry  using  explosive  charges  or 
"shots".  Departures  from  cylindrical  spreading  are  then  ascribed  to 
volume  attenuation  with  least  square  fits  made  to  the  loss  after 
removal  of  the  cylindrical  spreading  term. 

The  basic  assumption  of  this  technique  is  that  the  range  dependence 
of  the  spreading  losses  associated  with  axis-to-axis  transmission 
measured  using  bandpass  filtering  of  shot  signals  exhibits  a  1/R 
average  behavior.  This  assumption  has  been  justified  by  the  observa¬ 
tion  that  the  rms  sum  of  normal  modes  (representing  a  range  average 
of  a  CW  signal  level)  behaves  cyl indrical ly,  and  by  calculations  of 
DiNapoli  using  the  FFP  model  which  indicate  that  the  influence  of 
boundaries  on  axis-to-axis  propagation  is  felt  only  at  extremely  low 
frequencies. 

The  questions  left  unanswered  are: 

1.  What  frequency/range  dependent  effects  are  introduced  by 
the  shots  and  the  filtering? 

2.  What  frequency/range  dependent  effects  are  introduced  by 
range-dependent  changes  in  the  sound-speed  profile? 

While  geometries  for  experiments  are  chosen  with  great  care  to  avoid 
the  issues  raised  in  question  2,  data  from  the  KIWI  measurements  were 
shown  showing  significant  range  dependence  in  the  "absorption"  which 
appears  to  be  explainable  in  terms  of  unanticipated  intervening  bathymetry. 

It  has  also  been  pointed  out  by  Labianca  that  the  introduction  of  a 
complex  index  of  refraction  makes  the  wave  equation  non-self-adjoint 
and  appears  to  radically  alter  the  eigenvalue  spectrum.  (In  the  example 
of  a  surface  duct  the  spectrum  jumps  from  real -continuous  to  complex- 
discrete).  Whether  this  is  an  important  distinction  is  not  resolved 
but  the  inclusion  of  volume  absorption  by  a  complex  index  of  refraction 
may  cause  unexpected  problems. 


C.  Range  Dependences  in  the  Environment 

The  key  question  here  is  how  are  range  dependences  in  the  sound- 
speed  profile  and/or  bottom  profile  accomodated.  There  are  four  basic 
approaches  to  the  problem: 
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1.  The  parabolic  equation  technique  (BTL-1 )  which,  within  its 
small-angle  limitation  handles  range  variations  with  no  additional 
approximations. 

2.  The  NUSC-2  model  which  includes  only  first  order  mode¬ 
coupling  (i.e.  an  adiabatic  approximation)  but  also  includes  effects 
due  to  the  lack  of  cylindrical  symmetry, 

3.  Other  adiabatic  approaches  (NUC-1,  NRL-3)  where  only  the 
rms  sum  of  modes  is  considered. 

4.  A  non-adiabatic  technique  where  the  medium  changes  discon- 
tinuously  at  specific  ranges.  At  such  a  range  the  field  (as  computed 
from  the  modes  in  the  first  region)  is  expanded  in  terms  of  the  normal 
modes  of  the  new  region  (thus  giving  the  mode  amplitudes).  Fitzgerald 
of  NRL  is  using  this  technique  neglecting  all  backscatter  from  the 
discontinuity,  and  DiNapoli  is  considering  FFP  techniques  where  the 
backscattered  field  is  also  included  in  an  iterative  way. 

Only  the  parabolic  equation  contains  the  full  coupling  between  modes. 
Additional  work  at  CU  is  under  way  to  approximate  mode  coupling.  Holford 
presented  an  exact  solution  he  had  obtained  for  the  field  due  to  a  line 
source  (taken  to  be  the  z  axis)  for  a  biquadratic  variation  in  n2  with 
range  and  depth.  The  solution  is  written  in  terms  of  an  integral  of  the 
form 


*  *  g(C)eikf(x->,;E)  d5 


(20) 


While  this  solution  does  not  permit  boundaries  it  does  allow  for  a  range 
dependent  sound  channel  and  may  prove  quite  useful  as  a  control  case. 
This  result  is  documented  in  a  BTL  report. 


D.  Modeling  the  Ocean  Bottom 

The  basic  approaches  to  modeling  the  ocean  bottom  are  relatively 
straightforward  within  the  context  of  each  model.  Kutschale  and  DiNapoli 
have  devoted  considerable  attention  to  the  problem  of  layering  in  the 
bottom  and  its  implications  for  propagation.  Both  felt  that  contribu¬ 
tions  from  the  continuous  spectrum  in  modeling  the  response  from  the 
bottom  at  short  ranges  could  not  be  neglected. 

Some  of  the  normal-mode  programs  (NADC,  NUC-1)  allow  for  the  full 
computation  of  the  mode  amplitudes  for  a  layered  bottom  or  the  intro¬ 
duction  of  a  "reflection  coefficient"  as  a  function  of  angle  of  incidence 
for  plane  waves  at  the  angle  of  the  ray  equivalent  for  the  mode.  DiNapoli 
had  made  comparisons  of  FFP  with  the  CONGRATS  Ray-Tracing  Program  where 
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the  FFP  was  used  to  determine  the  plane-wave  reflection  coefficients 
for  a  multi-layered  bottom  which  were  then  used  by  CONGRATS  for  rays 
reflecting  from  the  ocean-bottom  Interface.  Comparisons  of  full  FFP 
and  CONGRATS  runs  showed  alarmingly  large  differences  in  the  bottom- 
bounce  regions.  This  also  raised  several  interesting  (and  unanswered) 
questions  concerning  the  reverse  problem  -  that  is  how  to  use  measured 
plane  wave  reflection  coefficients  (presumably  determined  with  the  aid 
of  ray-tracing)  in  wave  models. 


E.  Modeling  Source  Characteristics 

The  two  primary  questions  in  modeling  sources  are: 

1.  How  to  model  broadband  shots  (and  the  shot  processor)? 

2.  How  to  model  source  directionality? 

The  obvious  brute-force  way  to  model  broadband  sources  is  by  running 
normal -mode  programs  at  a  variety  of  frequencies  and  convolving  the 
result  with  the  source  spectrum.  DiNapoli  is  in  the  process  of  extending 
the  FFP  to  do  this  in  a  more  elegent  (if  not  more  efficient)  way. 

Source  patterns  can  be  modeled  with  multipole  expansions  (presently 
being  studies  at  CU),  though  this  involves  reworking  the  problem  for  each 
new  environment.  No  direct  weighting  of  mode  amplitudes  with  the 
directivity  pattern  for  the  ray  equivalent  has  been  attempted.  The 
parabolic  equation  technique  might  be  applied  to  this  problem  since  in 
its  transform  space  the  independent  variable  is  quite  similar  to  an 
angle. 


F.  Smoothing  Wave  Results 

The  question  was  posed  as  to  how  one  obtains  mean  transmission  loss 
with  normal  mode  programs.  DiNapoli  fits  the  results  of  the  FFP  to 
cyclindrical  spreading  plus  a  constant  which  is  determined  by  regression 
analysis.  This  technique  fails  if  there  are  strong  convergence  zones. 
Bartberger  has  an  extremely  simple  approach  which  is  to  use  a  weighted 
sliding  window.  Examples  were  presented  which  showed  that  as  the  width 
of  the  window  was  increased  the  rapid  oscftlations  were  removed  yielding 
results  very  similar  to  the  FACT  predictions. 


G.  Signal -Processing  Applications 

The  entire  field  was  not  considered  -  only  the  question  of  coherence 
and  wave-front  distortion.  Simulation  work  is  being  done  by  DiNapoli 
on  the  problem  of  wave-front  distortion  in  the  vertical  plane.  The 
problems  of  horizontal  coherence  across  arrays  are  considerable  and  at 
present  are  being  addressed  only  by  ray  acoustics.  Application  of  the 
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NUSC-2  Model  to  this  problem  might  be  fruitful  if  the  necessary  environ¬ 
mental  information  were  available. 


SUMMARY  AND  ACKNOWLEDGEMENTS 


From  the  point  of  view  of  AESD  the  workshop  was  quite  successful  in 
assessing  the  state-of-the-art  in  non-ray-tracing  acoustic  propagation 
models.  Subsequent  feed-back  from  participants  suggests  a  consensus 
view  that  it  was  successful  from  their  viewpoint  in  helping  both  to 
avoid  duplication  of  effort  and  to  illuminate  potential  areas  of  R&D 
rel event  to  the  Navy's  advanced  development  needs.  17  models  were  pre¬ 
sented  with  a  total  of  over  50  predictions.  The  success  of  this  work¬ 
shop  was  the  direct  result  of  the  enthusiastic  support  and  participa¬ 
tion  of  everyone  concerned. 
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